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Theory of Spontaneous Emission Noise in Open 
Resonators and its Application to Lasers 
and Optical Amplifiers 
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Abstract — A theory of spontaneous emission noise is presented based 
on classical electromagnetic theory. Unlike conventional theories of 
laser noise, this presentation is valid for open resonators. A local Lan- 
gevin force is added to the wave equation to account for spontaneous 
emission. A general expression is found relating the diffusion coeffi- 
cient of this force to the imaginary part of the dielectric function. The 
fields of lasers and amplifiers are found by solving the wave equation 
by the Green's function method. The lasing mode is a resonant state 
associated with a pole in Green's function. In this way, noise in lasers 
and amplifiers is treated by a unified approach that is valid for either 
gain guiding or index guiding. The Langevin rate equations for the 
laser are derived, the theory is illustrated with applications to trav- 
eling wave and Fabry-Perot amplifiers and Fabry-Perot lasers. Sev- 
eral new results are found: optical amplifier noise increases inversely 
with quantum efficiency; spontaneous emission into the lasing mode is 
enhanced in lasers with low facet reflectivities; and the linewidth of a 
Fabry-Perot laser with a passive section decreases as the square of the 
fraction of the cavity optical length that is active. 

I. Introduction 

THE TREATMENT of noise in semiconductor lasers 
has been quite successful in describing many proper- 
ties including the frequency spectrum arid amplitude dis- 
tribution of intensity noise, the linewidth and line shape 
of single-mode lasers and mode partition noise. Despite 
these successes, the current theories are based on assump- 
tions that are only approximately valid for most semicon- 
ductor lasers and amplifiers. The traditional theories of 
noise in lasers [l]-[3], begin by assuming that the laser 
can be treated as a closed cavity having a discrete set of 
orthogonal modes. The electromagnetic field is expanded 
in terms of these modes. In modeling the laser, facet losses 
are replaced by losses uniformly distributed throughout 
the cavity. This limitation of conventional laser noise the- 
ory is well known and has been addressed in a number of 
papers dealing with output coupling [4]-[6]. Neglect of 
facet losses is reasonable in the modeling of gas lasers 
with high reflecting facets, for which the theories of laser 
noise were originally developed. This approximation be- 
comes more questionable when it is applied to a conven- 
tional semiconductor laser having cleaved facets with only 
30-percent power reflectivities. 

The expansion of the electromagnetic field in terms of 
discrete resonator modes cannot be applied to traveling 
wave optical amplifiers, because, ideally, these devices 
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have no modes at all. In treating traveling wave and res- 
onant optical amplifiers, a different approach is used [7]- 
[9]. The electromagnetic energy in each transverse mode 
mode is considered to be locally excited and then ampli- 
fied during subsequent propagation in the device. 

Another well-known example of an open cavity, where 
conventional laser noise theory fails, is the gain-guided 
laser. This problem was solved by Petermann [10], who 
showed that gain guiding can greatly enhance sponta- 
neous into the fundamental transverse mode of a wave- 
guide. Streifer et al [11] applied Petermann's result to 
explain the iriultilongitudinal mode behavior of gain- 
guided lasers. However, Petermann did not derive his re- 
sult from a general theory of laser noise and for this rea- 
son he stimulated many other papers ; both defending and 
attacking him [12]. 1 A general treatment of traveling wave 
amplifier noise, incorporating Petermahn's result, has been 
given by Haus and Kawakami [13]. 

The purpose of this paper is to give a general formu- 
lation of spontaneous emission noise that is valid for open 
cavities and that can be readily applied to most semicon- 
ductor lasers and amplifiers. We will apply this theory to 
four cases: a traveling wave amplifier, a Fabry-Perot am- 
plifier, a Fabry-Perot laser with arbitrary facet reflectiv- 
ities, and the calculation of the Lorentz linewidth of a 
Fabry-Perot laser with a passive section. The theory is 
presented in a form suitable for application of other de- 
vices of current interest such as distributed feedback and 
coupled cavity lasers. 

Our starting point is the same as that of Landau and 
Lifshitz, who treat electromagnetic fluctuations in dielec- 
tric media [14j. We assume that a field source (Langevin 
force) FJX) accounts for spontaneous emission by the 
carriers at angular frequency o>. We find the mean square 
(diffusion coefficient) of the Langevin force, which is pro- 
portional to the imaginary part of the dielectric function 
e w (x). 

When this source of spontaneous emission is included, 
the wave equation becomes inhomogeneous. The fields 
generated by this source may be calculated by the Green's 
function method. The Green's function depends on device 

'M. Newstein. in [12], has criticized Petermann [10]. His criticism has 
to do with how spontaneous emission into the lasing mode is denned. The 
definition of spontaneous emission rate R used here in (48) leads to no 
disagreement with Newstein. 



0733-8724/86/O3O0-0288$01.O0 © 1986 IEEE 



HENRY: EMISSION NOISE IN OPEN RESONATORS 



289 



geometry and relates the universal noise source to the perr 
formance of specific devices. 

Our treatment of laser noise is semiclassicai. The quan- 
tum theory of radiation is hot used and the amplitude of 
the noise source is determined merely by requiring that 
the spontaneous emission generated be exactly that re- 
quired by statistical mechanics in the idealized case of an 
optical waveguide in equilibrium with a semiconductor. 
This approach is close to one we used earlier to relate the 
coefficients of absorption, gain, and spontaneous emis- 
sion [15]. 

These methods are well known in electromagnetic the- 
ory, but have not been used very much in the description 
of lasers. Petermann's [7] calculation is essentially the 
calculation of a Green's function for a uniform wave- 
guide. Our source of spontaneous emission resembles that 
of Haus and Kawakami [13], but we do not use an effec- 
tive temperature as they do. 

H. Spqntaneous Emission Noise 

A. Basic Equations 

Our starting point is the wave equation for the electro- 
magnetic field. For simplicity, we will avoid the use of 
vector fields. We assume that in the laser or amplifier, 
waves propagate primarily along the z axis. E(x, t) is the 
transverse electric field component with polarization vec- 
tor in the x-y plane. If we describe the real electric field 
in terms of its complex Fourier components EJx) 



E(x 



, t) = \ E„(: 
Jo 



x) exp (-iur) du + c.c. (1) 



where c.c. is the complex conjugate, the wave equation 
for EJx) is 

V 2 E„(x) +-jDJx) = 0. (2) 

The electric displacement vector D w is assumed to be 

D w (x) = e^x) + K„(x) (3) 

where the first term is due to the displacement induced by 
the electric field and the second term is a Langevin force 
associated with spontaneous polarizations set up by the 
semiconductor carriers [14]. The average value of K„ is 
zero. 

Substitution of (3) into (2) gives 

,.2 



(4) 



(5) 



where Langevin force F w (x) is defined as 
FJx) m —jKJx). 

B. Green 's Function Solution of the Inhomogeneous 
Wave Equation 

Equation (4) may be solved by first finding the Green's 
function G w (x, x' ) defined by 

(v 2 + J e w ) Gjx, x') = 6(x - x') (6) 



where 5(x - x' ) is the Dirac delta function. Then the field 
generated by spontaneous emission (and amplification) is 
given by 



E a (x) = j dx' G w (x,x')F w (x'). 



(7) 



For brevity, we use x = (x, z), where x represents both 
lateral coordinates x and y. We will assume that there is 
a complete set of transverse modes 0„(x) satisfying 

r 22 '2 



^2 + f *„(x) = fy H (x). 

The completeness relation is 
v^(x) <t>„(x') 

n 

where 



= 5(x-x') 



(8) 



(9) 



j <t> n <t> m dx m <<M> m > = ««<*A>. 
In the case of index guiding, <j>„(x) is real and <<M>«> = 

Making use of the completeness relation, we can ex- 
pand the Green's function in terms of 4> n 



G„(x, X') = Zgn(Z 9 Z') 



<t> n {x') 



(10) 



The one dimensional Green's function g„{z> z' ) satisfies 
g„{z, z') - 6(z - z') ' (11) 



dz 2 + K 



as can be seen by substituting (9) and (10) into (6) and 
using (8). In some applications, such as the distributed 
feedback laser, k n will also depend on z. The solution of 
(11) has been discussed at length by Morse and Feshback 
[16]. It is given by a one-dimensional Green's function 



gn{z> z') 



Z w+ fe>) Z„_(z<) 
W m 



(12) 



where z> and z< are the greater and lesser values of z and 
z' and Z n +(z) and Z n .(z) are the solutions of the homo- 
geneous part of (11), satisfying the boundary conditions 
for positive and negative z, respectively. The Wronskian 
W n is defined by 



(13) 



where Z' n = dZJdz. It can be shown [15] that W n is in- 
dependent of z. the function g„(z, z') has a slope discon- 
tinuity allowing it to satisfy (11): 

C. Laser Amplifier 

We will be primarily interested in fundamental single 
transverse mode operation of amplifiers and lasers. In this 
case, we need only retain the n = 0 term in (10). For an 
amplifier extending from z = 0 to z^= L, the field at L is 
given by 
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EJx, L) = Zo + (L) *o(x) 



(14) 



where 



{Zo^oFJ = J Zo- A j <ir F w (x, *). 
ZX Layer 

In the case of a laser operating in the fundamental trans- 
verse mode, it is useful to think of W 0 as a function of the 
real angular frequency o>, which can be analytically con- 
tinued to complex values g>. The zero's of W 0 (w) are the 
complex values £>, . These values uy { are the poles of EJx) 
in (14), which correspond to the longitudinal modes of 
the system. We see from (13) that at ^(aj,) = 0, the func- 
tions 2q + (z) and Zq_(z) are proportional to one another. 
We will limit our discussion to single-mode lasing oper- 
ation, i.e., we will assume that there is only a single pole 
o) 0 near the real axis. We can define the functions Zo+ and 
Zo- to be equal at this pole: Zq + (z) = Zq-(z) = Zq(z) for 
w near £> 0 . For a laser near threshold and for frequencies 
near a> 0 , we can expand W 0 as W 0 (w) = (dW 0 /do)) (a? - 
d> 0 ). The field (14) then becomes 



EJx) = 



Zpfe)<fro(*) <Zo0 o F„) 
<<M>o> ~!~r (w - a> 0 ) 



(15) 



As the carrier number in the laser increases and thresh- 
old is approached both the real and imaginary parts of a> 0 
change. This is illustrated in Fig. 1. The pole wq can be 
expressed as 

aAG iAG iAG 
wo = o>o + -y- + -y- = wq + — (1 - ice) (16) 

where o> 0 is the cavity resonance at threshold, AG is the 
net gain (s~ l ) and a is the linewidth parameter [17], which 
determines the rate of change of the cavity resonance with 
gain. The gain AG determines the rate of growth or decay 
of the field intensity in the absence of spontaneous emis- 
sion. Below threshold, AG is negative. In a laser that is 
uniform along z, oc is the ratio of the changes in the real 
and imaginary parts of the refractive index with carrier 
number [17]. Henry and Kazarinov [18] have used com- 
plex resonance frequency (16) to discuss mode suppres- 
sion in cleaved coupled cavity lasers. 

In a laser below threshold, AG, which is related to the 
carrier number N, can be taken to be a fixed parameter. 
Fluctuations in N can be neglected. Above threshold, TV 
is affected by light intensity and the laser field is very 
sensitive to changes in N. All phenomena that distinguish 
the laser oscillator from a resonant amplifier, e.g., gain 
saturation (Fermi-level pinning); the change in the light 
intensity distribution from exponential to Gaussian; relax- 
ation oscillations; and the enhancement of laser iinewidth 
are due to. the dynamic interaction of the carrier number 
and the light intensity. This interaction is most conve- 
niently described by a set of coupled differential equations 
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Fig. 1. The pole associated with the lasing mode in the complex w plane. 

for field amplitude 0(f) and carrier number N. The com- 
plex equation for 0(f) is often transformed into two real 
equations for light intensity /(r) and phase 0(f). 

We will now derive the equation for the complex field 
amplitude 0(f). We note that the spatial dependence of 
EJx) in (15) is given by Zo(z) 4> 0 {x) and define a field am- 
plitude 0 W as 

K(x) = B^Zo(z) Mx). (17) 

Using (1) and (17), we can write the time dependent elec- 
tric field £(x, r) of the laser as 

E(x, f) = 5/5(0 Zoiz) <*>o(*) + c.c. (18) 

where the field amplitude 0(f) is given by 



0(f) = \ 0 W exp (-iwf) dw. 
Jo 



(19) 



The constant B in (17) is for normalization. We choose 
its value (in Section V-B) so that the instantaneous num- 
ber of photons in the cavity /(f) (integrated electromag- 
netic energy divided by fiw 0 ) is given by 

/« = 0(0 0(f)*. (20) 

An equation for 0(f) can be determined by multiplying 
(15) by -i fo (co - a> 0 ) exp (-/«/) dw and using (16), 
(17), and (19). The result is 



0 - J-/o> 0 + (1 - id)] 0 + F^t) 



where 



and 



-iiZrfoFjx, Q) 

£<<Mo> ;p 
do) 



F(x, f) = f FJx) exp {-iut) da. 
Jo 



(21) 



(22) 



(23) 
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ACTIVE 
L — 



PASSIVE 
L, 



Fig. 2. Diagram defining the Fabry-Perot cavities considered in this pa- 
per, (a) uniform cavity, (b) cavity with passive section. 

E. Green 's Function for an Infinite Waveguide 

We will now give several examples of Z n+ |,Z„_ and W n 
that are useful later in this paper. For an unending wave- 
guide (or a traveling wave amplifier with perfect antire- 
flection AR coatings), we can choose 

Z n± = exp(±i*„z) (24) 

then the Wronskian (13) is 

W n = 2ik n (25) 

and the one-dimensional Green's function (12) is 

(26 ) 

The complete Green's function (10) is 

Lr^X, X ) - Zj . (27) 

n 2ik n {<t> n <t> n } 

F. Green '$ Function for the Fabry-Perot Cavity 
Consider a cavity extending from z = 0 to z = L with 

field reflectivities r t and z - 0 and r 2 atz = L (Fig. 2(a)). 
It is easily shown that for the fundamental transverse mode 

Zq_ = r, exp (ikoz) + exp (-&oz) (28) 

satisfies the boundary condition at z = 0 and 

Zq + = r,[exp (ik 0 z) + r 2 exp (2i* 0 L - i*bz)] (29) 

satisfies the boundary condition at z = L. 

The factor r, is added to (29) so that at threshold, when 
r x r 2 exp {2ik^L) =.1, Zo+ 2q_. The Wronskian (13) 
for the fundamental mode is given by 

W Q = 2j*br,[l - r x r 2 exp (lik^L)]. ' (30) 

Near threshold *o = Jfc 0t h + A&o and . 

1 - r,r 2 exp (2i*oL) 



1 - exp (2/A*oZ,) ~ -2/A^L (31) 



where 



(32) 



and is the group velocity. Substitution of (32) with AG 
■ *>*o Ago into (30) yields 



0) - co 0 



4* 0 L/' 1 

Therefore, for the Fabry-Perot cavity 

dW 0 = 4k 0 Lr [ 
aw 



AG ] 



(33) 



(34) 



HI. Diffusion Coefficients of the Langevin 
Forces 

A. Diffusion Coefficients of F u (x) 

it can be generally established [19, sec. 118] from the 
requirement of stationarity, (i.e., that <F(x, t) F(x\ t + 
t)> is independent of r), that 

<W F w ,(x')> - <F w (x)» F tt <x')*> - 0 (35) 

and that < F w (x) F w <*')*> is real and only nonzero when 
ci>' = «. We assume that the Langevin forces have neg- 
ligible spatial correlation. (They should only be corre- 
lated over the distance a carrier can travel before being 
scattered, * 10~ 2 pirn). Therefore 

<F w (x) F w ,(x')*> = 2D FF *(x) 6(x - x') 5(o> - o>') (36) 

where the 4 diffusion coefficient" D FF *(x) is real. 

We can determine D^x) by requiring light traveling 
in a waveguide many absorption lengths long and having 
only bandgap absorption losses comes to equilibrium with 
the semiconductor. In equilibrium, the average noise 
power P N in the z direction, within angular frequency in- 
terval Aw will be L n v gn hu < * w > multiplied by the number 
of modes per unit length (dkjdw) (Ao>/2?r), hence 



o/ = " £j 

2T n 



(37) 



where v gn = dw/dk n is the group velocity of the nth mode, 
Afc„ is the change in wave vector for change Aa>, and the 
sum is over the transverse modes. In equilibrium, the 
chemical potential of the photons will be eV, the separa- 
tion of the quasi-levels of the conduction band and val- 
ence band of the semiconductor. Then the change in free 
energy during emission or absorption is zero [15], With 
this chemical potential, the average mode occupation 
number {19, sec. 53] is 



exp 



(38) 



1 



In equilibrium, hu> > eVand .(n„) is positive. In a laser 
h<j> 0 < eV, < nj> is negative and the positive quantity n sp 
is normally used [16]. 

We will now derive an expression for P N (37) by use of 
the Green's function. Comparison of the two expressions 
for P N will allow us tp determine the diffusion coefficient 
D FF *(x) (36). The power P N is given by an integral over 
the z component of the Poynting vector S = cE x H/4tc 
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• exp (-/(oj - a>') r) + cc.)). (39) 
From Maxwell's equation 

-ic d£L 



a) dz 



(40) 



We will assume that our ideal waveguide is index guided 
so that <t>„ is real and <<M>n> = 1 ■ Writing E w as a Green's 
function (7), (10), (12), and using (36), we find 



£. Diffusion Coefficient for F & (t) 

The Langevin force Fp(r) acting as a source of sponta- 
neous emission into the lasing mode at co 0 is defined by 
(22) and (23). Using these equations, (F fi (t) F^{t')) is 
given by 

W)F^')*> = 

( du exp (-ia>(r - r') (Zq^D^ZJ^) 
Jo 



(£Z„ + (L) Z„ + (L)'* + c.c.) f \Z„.\ 2 
Jo 



j2 <fc. (41) 

Evaluating Z n± and W n with (24) and (25), we find 



" 16™ n fcfc? 



(42) 



where k n = k' n + ik„. We have neglected k; 2 compared to 
k' n 2 and approximated 1 - exp (-2/c^L) ~ 1. The later 
approximation is valid for a waveguide assumed to be 
many absorption lengths long. Using the equation for the 
lateral field (8) it is readily shown that 



2KK 



(43) 



where e" is the imaginary part of the dielectric constant 
e. Therefore 



P - A ^ 4 V <0„2ZW n > 
N 2* 4o> 3 n <<t> n e"<t> n > ' 



(44) 



This can be made to agree with the statistical mechanical 
result (37) if 

2Dfr{x) = jr^in^y (45) 

A corresponding expression was derived by Landau and 
Lifshitz using quantum mechanics for thermal fluctua- 
tions of the electromagnetic field [14, eq. (88.11)]. Our 
expressions agree if we set the bias voltage V = 0 and 
add a term for the zero-point oscillations of the field by 
changing <n w > to (nj) + \. Equation (45) is a form of 
the fluctuation dissipation theorem [18, sec. 123]. It ex- 
presses the fact that in equilibrium dissipation by optical 
absorption must be balanced by spontaneous emission 
causing field fluctuations. Having determined 2D FF * (x), 
we can then apply it to lasers and amplifiers which are 
generally not in equilibrium. Expressing (co/c) e"(*) ~ 
-*'(*) gW, where n' is the real part of the refractive 
index and g(x) is the gain, (45) can be written as 

2D^(x) = ^ n'(x) g{x) n sp (46) 
where n sp is the spontaneous emission factor (38). 



B 2 i<<M>o>| 2 gf| 



(47) 



where the angular bracket on the right side is an integral 
over jc and z. We are only interested in the frequency com- 
ponents of oj near u> 0 . Therefore, we will approximate 
IDff* by its value at a> 0 and remove it from the integral. 
The remaining integral is 27r6(r - t f ). This approxima- 
tion only amounts to replacing a narrow function by a delta 
function. Therefore 



(48) 



where 



a = — — p . (49) 

We will show later that R is the average rate of sponta- 
neous emission into the mode. A similar calculation, 
making use of (35), shows that 

<Ff£t) Fe(t')> = <F^(r)* F,(f')*> - 0. (50) 

To complete the evaluation of R, the normalization con- 
stant B must be chosen. We choose B so that when the 
field is expressed in terms of 0, the number of photons in 
the cavity is / = |£| 2 . With this choise, the field £ is the 
classical average of the photon annihilation operator b 
[20]. This choice is convenient because it allows the Lan- 
gevin equations for light intensity / and carrier number N 
and the related diffusion coefficients to be expressed sim- 
ply in terms of parameter R. In a dispersive medium, the 
electromagnetic energy per unit volume U is given by [14, 
eq. (61.10)] 



8* 



dcc 



(e'a>) E 2 + 



4 



Using (18) and (40), this becomes 



U = 



4?r 



f- 



(o>/«' 2 )|2y 2 + -|Z6| 



(51) 



|/3| 2 - (52) 



At any point in the cavity, Z consists essentially of back- 
ward and forward propagating waves with propagation 
constant u>n'lc. In this case 



(53) 
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where only oscillatory terms, which contribute negligibly 
to the integral over the cavity, length, have been ne- 
glected. Combining (52) and (53), we find 



U 



B 2 <fott\Zo\ 2 nfr t 
2tt 



I/3I 2 



(54) 



where n g = w(dn'/da)) + ri is the group index. Equating 
$ Udx/huQ to |(3| 2 , we find 



B 2 = 



2ttHojq 



<4>Hq> <Zff/i6%,Zo> 



(55) 



where the bracket indicates an integral and (0o*'"*4ff) 
= flo*g<<M>(f>- Substitution in (49) results in 

4K<4<Z$gonon sp Zo) {Z$n^ n^2^> 

(55) 



R = 



where £ is the gain guiding enhancement factor of Peter- 
mann [10] 



K = 



<Wo) 2 
|<*o*o>l 2 



(57) 



IV. Application to Amplifiers 



2ir(g 0 ~ oto) 
The power amplification of the signal is 
A = exp [(g 0 - oq) L]. 



#0 



2th? 



■04 - 1). 



(61) 



rivations use ij = 1, thereby underestimating the noise 
figure by about x2, 

B. Fabry-Perot Amplifier . 

In practice, end reflections are unavoidable and every 
TWA is actually a Fabry-Perot amplifier (FPA). We will 
calculate the noise power of an FPA under the same as- 
sumptions as made for the TWA: fundamental transverse 
mode operation, an additional distributed loss a 0 and in- 
dex guiding. The noise power is calculated exactly as in 
the case of a TWA, except that Z 0± and W 0 for the Fabry- 
Perot cavity (28)-(30) are inserted into (41). For simplic- 
ity we will regard the two facet field reflectivities as real 
and positive. The only effect of complex reflectivities is 
to introduce phase shifts which slightly alter the effective 
cavity length L. The integral of |Zo- 1 2 from 0 to L is given 
by 

<|Zo_| 2 > = [r 2 + exp (-(go - *o) D\ 



A. Traveling Wave Amplifier (TWA) 

We will calculate the noise power P N generated by 
spontaneous emission in a TWA having perfect antire- 
flection facet coatings, assuming propagation in the fun- 
damental mode and index guiding. The steps of calcula- 
tion are the same as those used for deriving (44) with only 
one transverse mode n = 0 included. An additional loss 
per unit length ot 0 is assumed, so that kg = (- 1/2) (g 0 - 
a 0 ). The diffusion coefficient Dpf*(x) in (44) is given by 
(46). The resulting noise power P N is given by 



[exp ((g 0 - op) L) - I]. (58) 



(59) 



In a laser having no leakage currents, the quantum effi- 
ciency is given by 



(60) 



It is the ratio of stimulated emission not dissipated to total 
stimulated emission. With these definitions we can write 



We see that two nonideality factors increase noise: n sp and 
7j. In good, long wavelength lasers n sp » 2 and tj » 0.5, 
so that these factors increase the noise by about x4. 

Equation (61) is well known (See [7]-[9], [13] and [21, 
references].) However, as far as we know, previous de- 



exp [(g 0 - a 0 ) L] - 1 
go ~ <*o 



(62) 



and 



iZo+(L) Zo + (L)'* + ex. ~ r?(l - r\) 

' exp [(g 0 - oo) L] 
Upon substitution into (41), we find 



(63) 



sp 



[exp ((g 0 - oo) L) - 1] 



2ttj 

f[rj exp ((g 0 — a) L) +.!](! - rjfl 
I |1 - r,r 2 exp (2^oL)| 2 J 



(64) 



This differs from P N for the TWA (58) by an additional 
resonant term in brackets which goes to 1 as r, and r 2 go 
to 0. With tj = 1, (64) agrees with the expression given 
by Yamamoto [22] based on the earlier work of Gordon 
[8]. It is easily shown that optical power is amplified in 
passing through a FPA by 



A z= 0 ~ r 2 ) (1 - rf) exp [(g 0 - op) L) 



|1 - r,r 2 exp (2i^L)| 2 



(65) 



Therefore, 



[1 - exp (-(g 0 - oo)L)] 



2T7J 

[r] exp ((gp — otp) L) + 1)] 



r- r\ 



Near threshold r x r 2 exp [(g 0 - a) L] ** 1 and 



P N 



r,r 2 ) (r, + r 2 ) 

A ' 



2nj(l - rf) r 2 



(66) 



(67) 



C. Enhancement of Amplifier Noise by Gain Guiding 

In the case of gain-guided fundamental-mode opera- 
tion, we obtain an increase in noise power flux P N by a 
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factor K (57) which is exactly the result obtained by Pe- 
termann [10]. This is easily seen by examining (14) for 
the field of the fundamental mode at the end of an ampli- 
fier. Calculation, of the noise power from (14) and (39) 
involves calculating 



N |<4o*o>| 2 



- "oSo (68) 



Ptf(OUt) = 



(69) 



The above analysis is much simpler than that of Ya- 
mamoto [22] , who calculated noise of an amplifier chain 
from the viewpoint that the optical signals and noise are 
collections of photons having shot noise. In this paper, 
we emphasize the wave viewpoint. An analysis of an am- 
plifier chain, similar to that given here has been made by 
Henry [24]. 

V. Application to Lasers 

A. Langevin Equations for land <t> 

The Langevin force equation for a single mode laser 
was derived earlier (21). The nonzero diffusion coefficient 



of the Langevin force F 0 (t) is determined by a single real 
parameter R given by (56). The fields are often more con- 
veniently expressed in terms of intensity / = : j30* and 
phase </> 

P = l m exp [-/(a> 0 f + </>)] 



or 



|<<M>o>! 2 "' 

where the last equality is easily derived from (8). Thus 
P N is proportional to K which goes to 1 in the case of 
index guiding. 

D. The Distribution of E u Emerging From a Chain of 
Optical Amplifiers 

In calculating error rates in optical communications, the 
statistical distribution of the optical field at the detector is 
needed. The field (14) associated with amplified sponta- 
neous emission E w is the sum (an integral) of many con- 
tribution F w (x). According to the central limit theorem 
[23] the distribution of the sum of many statistically in- 
dependent random variables with common distributions is 
a Gaussian. Therefore, the real and imaginary parts of E u 
have Gaussian distributions. By applying (35) it is easily 
seen that the real and imaginary parts of E u have identical 
(statically independent) Gaussian distributions. That is, if 
E w can be represented as a vector in the complex plane. 
The distribution of this vector is isotropic and along any 
direction (fixed phase angle) the distribution is the same 
Gaussian. 

When passed through subsequent linear amplifiers 
(having no gain saturation) the distribution of E u will be 
unchanged except in magnitude, but each amplifier in the 
chain will generate additional fields. It is well known that 
the sum of several Gaussian random variables is also a 
Gaussian random variable. Therefore, E™ at the end of 
the amplifier chain will also have an isotropic Gaussian 
distribution. The power reaching the detector will be pro- 
portional to < | E™ 1 \ 2 ) arid of course it will be just the sum 
of amplified powers generated by each amplifier. It is eas- 
ily shown that for a chain of M amplifiers of amplification 
A, in which each amplification is followed by an atten- 
uation A~. x so that the net amplification of the chain is 
unity that the noise power at the output is 



<t> = -wo' + - On 0 - In 0*). 



(70) 



(71) 



The Langevin equations and associated diffusion coef- 
ficients are completely determined by the transformation 
equations (20) and (71) which define / and <j> in terms of 
0, 0*, and t. Lax has given convenient and general for- 
mulas for obtaining the transformed Langevin equations 
and diffusion coefficients [25]. Making this transforma- 
tion, we find 



0 = 


i aAG + F4f) 


(72) 


/ = 


AG/ + R + F,(f) 


(73) 


2D„ = 


2RI 


(74) 


2D^ = 


R 

21 


(75) 


2D» = 


0. 


(76) 



Equation (73) shows that the average rate of spontaneous 
emission into the mode is R. 

B. Power Flow in a Laser 

The power P flowing along the z direction in a laser is 
given by integrating the Poynting vector over x. Express- 
ing E with (18), if with (40) and dropping 0 2 and 0* 2 
which oscillate at twice the optical frequency, we find 



SB 2 

P = — (&0Z6* + c.c.) /(r) <«o#>. 



(77) 



Eliminating B 2 with (55), we have 

2<zjr«6«*>Zo> () - (78) 

In an open cavity, the definition of 7 is ambiguous be- 
cause it depends upon what we choose to be the cavity 
length. The minimum extent of the cavity would have to 
include all active regions. If we choose to include passive 
regions beyond this length, 7 and <Z 0 /i6rt /r oZJ> will in- 
crease, however, physical quantities such as P involve the 
ratio of these two quantities and are unambiguous. 

C. Fabry-Perot Cavity Laser 

In lasers such as the Fabry-Perot cavity laser g 0 n 0y and 
VgQ are constant along the cavity length and R (56) reduces 
to 



R = 



4^0/WlZbl 2 ) 2 



do) 



(79) 
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0.1 02 0.3 0.4 

Fig. 3. Spontaneous emission rate R as a function of facet reflectivity. 

where ko = noo) 0 /c. We can evaluate R for the Fabry - 
Perot laser with the aid of (34) and (62), we find 

where and r 2 are the facet field reflectivities (Fig, 2(a)) 
and G 0 = £o y go- 1° th e case of equal facet reflectivities r { 
= r 2 = R m 

* " r L ( B(i5u)J Gon *°- (81) 

For index-guiding (^ = 1) and high-reflecting facets (r, 
— r 2 = 1)j R — QoHspOi which is the conventional result 
[17], [26], For R m = 0.3, corresponding to uncoated 
cleaved facets, the correction (bracketed term in (81)) is 
1.13. For R m = 0.01, the correction is 4.62. Hence the 
correction is modest for a conventional laser with un- 
coated cleaved facets and becomes quite significant for 
low-reflecting facets. The correction factor is plotted in 
Fig. 3. 2 

The facet powers of the Fabry-Perot laser can be easily 
found by evaluating P (78) at z = 0 and z - L with the 
substitution of (28) and (62). At z = 0, P x = -P is 

and at z = L\ P 2 = P is 

&-^ L^a- J (83) 

where oc m s L" 1 In ((rjr 2 ) 1 ). When r, = r 2 , the brack- 
eted terms reduce, to \. 

D. Linewidth of Fabry-Perot Laser Containing a 
Passive Section 

The broad linewidth of semiconductor lasers can be re- 
duced considerably by adding a long passive section to 
the laser cavity. We will calculate the linewidth reduction 
for an ideal Fabry-Perot cavity with a passive section 
shown in Fig. 2(b) in which there is a nonreflecting tran- 
sition between the active and passive sections. 

In this case, the Wronskian for the fundamental trans- 

2 This same correction factor was found by Ujihara 127) using a different 
approach. 



verse mode is 

W Q = 2/ViU - nr 2 exp (2i*oLo + 2ifc,L,)]. (84) 
It agrees with (30) except for the additional term 2i&,L ] . 
At threshold <o - w 0 and r,r 2 exp (HI^Lq + 2 /£,£,) = 1 . 
We can expand W 0 about a> 0 as before, see (30)-(34), but 
including an additional term 2Afe,L, = 2(w - w 0 ) Lj/i^, 
associated with the passive region. Expanding W 0 (w) as 
dW 0 /d& (o) - o) Q - /AG(1 - iot)l2) y we find that com- 
pared with the values for a uniform Fabry-Perot cavity 
(dW 0 /doj) increases by £ and AG decreases by £ , where 

t = f85 x 

* UJv# + L\fv gX { ** } 

is the fraction of the cavity optical length occupied by the 
active region. 

The function Zq(z) (28) is unchanged by the addition of 
the passive section. If we define / as the number of pho- 
tons in the active section, the spontaneous emission rate 
R (79) is unchanged except for a decrease by £ 2 due to the 
change in dW 0 Idw. Furthermore, with this definition of /, 
the relation of / and output power P x (82) is unchanged. 
The Lorentzian linewidth of the laser Av is given by the 
(spontaneous emission rate) (1 + a 2 )/47r/ [17], [26]. In 
this case, we have 

R(\ + a 2 )£ 2 

(86) 



Av = 



4tt1 



where R is given by (80), (81). Using (82) to relate / to 
the power out of the active end P, , we have 

V 16irLoP, 



| >i ± r z )(\ - r,r 2 )(l - r]) 
L r]r 2 In in 1 r 2 l ) 



(87) 



(1 + a 2 ). 



For equal reflectivities n = r 2 = R]? 9 the linewidth re- 
duces to 



(1 - *J 2 



R m In (I?- 1 )] 



(1 + a 2 ). (88) 



Compared with the Lorentzian linewidth calculated in 
[17], [26] for a uniform laser, Ay is increased by the en- 
hancement in spontaneous emission (81) and decreased by 
£ 2 . Line narrowing by a passive section can be substan- 
tial. For Lq = 300 fim, v g0 = c/4, L t = 3 cm, and v gl - 
cl.5, we find £ = ^ and £ 2 = t*W. This value of £ will 
reduce a linewidth of 100 MHz to 67 kHz. 

Summary and Discussion 
A theory of spontaneous emission noise in semiconduc- 
tor lasers and amplifiers is presented. The presentation is 
self-contained and based on classical electromagnetic the- 
ory. A local source of spontaneous emission in the form 
of a Langevin force F u (x) is added to the wave equation. 
The diffusion coefficient for this noise source is found by 
requiring consistency with the results of equilibrium sta- 
tistical mechanics. It is shown to be proportional to the 
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imaginary part of the dielectric function e£(x) and the 
spontaneous emission factor n sp (38), determining the de- 
gree of inversion of the amplifying medium. 
The Green's function method is used to relate the field 
t0 the noise source. In this way a theory of spon- 
taneous emission noise is developed that is valid for open 
resonators. It differs from conventional laser noise theory 
which relies upon expansion of the electromagnetic field 
in terms of the normal modes of a closed resonator [1]- 
[3]. 

The theory is applied to amplifiers and expressions are 
derived for the power spontaneously emitted from trav- 
eling wave and Fabry-Perot amplifiers. These expres- 
sions are conventional except that they contain an excess 
noise factor ij" 1 , where 17 is the quantum efficiency of the 
amplifier due to nonbandgap absorption and scattering 
losses. This factor appears to have been missed in earlier 
treatments of optical amplifiers. It is shown that the dis- 
tributions real and imaginary parts of the field generated 
by amplified spontaneous emission are identical Gauss- 
ians and that this result is valid for a single amplifier or a 
chain of linear amplifiers and attenuators. 

the lasing mode is shown to be associated with the pole 
of the Green's function. In this way, the same methods 
can be applied to treat lasers and amplifiers. The Lan- 
gevin rate equations for the laser field, intensity and phase 
are derived. The various diffusion coefficients of the Lan- 
gevin forces are shown to depend on a single parameter 
/?, the spontaneous emission rate. A general formula for 
R is given (56). This formula is evaluated for the Fabry- 
Perot cavity. It is shown that for a Fabry-Perot cavity 
with high-reflecting ends, R reduces to the conventional 
value Gpn S p, where G 0 is the gain per second. For lower 
values of reflectivity corresponding to an open cavity, R 
is enhanced (Fig. 3). 3 A Fabry-Perot cavity laser with a 
passive section is also treated. It is shown that the line- 
width of such a laser can be reduced by £ 2 where £ is 
given by the fraction of the optical length of the cavity 
occupied by the active region (85). 

Our theory is valid for gain-guided as well as for index- 
guided lasers. We find that for gain-guided single-trans^ 
verse mode lasers and amplifiers, spontaneous emission 
is enhanced by the K factor of Petermann [10]. 

In conclusion, we have presented a self-contained the- 
ory of laser noise, valid for open resonators, that unites 
the treatments of amplifiers and gain guiding with the 
conventional treatment of lasers. 

Nomenclature 
£(x, 0 Electric field. 

Eu(x) Fourier component of elec- 

tric field. 

D a {x) Fourier component of elec- 

tric displacement. 

K u (x) Contribution to Z> w by spon- 

taneous polarizations. 

'(After completing this paper it was found that Ujihara [27] had found 
the same effect.) 



F w (x) = (-o^/c 2 ) K a (x) 

e w (r)e'(jc) + fc'(r) 
C w (x,jc') 

g n (z, z') 
Z n± (z) 

AC 

Of 

Zo(z) 

ft, 
B 

F e (t) 

D FF »(x) 
(O = ~n tp 
Pn 

n'(x) + in"(x) 

g(x) = (-o)/c) n(xy 

rto» 8o> a o 

G 0 = goVgo 
R 

A 
L 



Langevin force of the wave 

equation (4). 
Dielectric function. 
Greens functions satisfying 

(6). 

nth transverse mode. 
Propagation vector of nth 

transverse mode (8). 
Green's function for the z 

coordinate of mode n 

(11). 

Solutions of homogeneous 
part of (11). 

Wronskian (12), (13). 

Pole associated with the las- 
ing mode (16). 

Threshold angular fre- 
quency of fundamental 
mode (16) . 

Net gain (s~') of fundamen- 
tal mode (16). 

Linewidth parameter of fun- 
damental mode (16). 

z dependence of fundamen- 
tal lasing mode (17), (18). 

Amplitude of the laser mode 

(18) . 

Fourier component of &(t) 

(19) . 

Normalization constant (17). 
Intensity of lasing mode 

(20) . 

Langevin force driving laser 
field amplitude (23), (24). 

Derivative of W Q used to ex- 
pand W 0 near laser mode 
(15). 

Group velocity of the fun- 
damental mode. 

Diffusion coefficient of 
F„(x) (36). 

Equilibrium mode occupa- 
tion number (38). 

Noise power of amplifier 
(42). 

Refractive index. 

Gain per centimeter. 

Refractive index, gain, and 
loss of the fundamental 
mode. 

Gains/s of fundamental 
mode. 

Spontaneous emission rate 
(56). 

Power amplification. 
Length of Fabry-Perot cav- 
ity (Fig. 2(a)). 
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K 



<t> 



U 



Field reflectivities of Fa- 
bry-Perot cavity (Fig. 2). 

Facet power reflectivity for 
n = r 2 . 

Enhancement of sponta- 
neous emission by gain 
guiding (57). 

Quantum efficiency (60). 

Phase of laser field (70). 

Diffusion coefficients of in- 
tensity and phase for laser 
(74)-(76). 

Electromagnetic energy 
density (51), (54). 

Fraction of optical length 
occupied by active region 
((85), Fig. 2(b)). 

Active and passive lengths 
in Fig. 2(b). 
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